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Abstract 

We report neutron diffraction and magnetization studies of the magnetic order in multiferroic 
BiFeOa- In ferroelectric monodomain single crystals, there are three magnetic cycloidal domains 
with propagation vectors equivalent by crystallographic symmetry. The cycloid period slowly grows 
with increasing temperature. The magnetic domain populations do not change with temperature 
except in the close vicinity of the Neel temperature, at which, in addition, a small jump in magneti- 
zation is observed. No evidence for the spin-reorientation transitions proposed in previous Raman 
and dielectric studies is found. The magnetic cycloid is slightly anharmonic for T=5 K. The an- 
harmonicity is much smaller than previously reported in NMR studies. At room temperature, a 
circular cycloid is observed, within errors. We argue that the observed anharmonicity provides 
important clues for understanding electromagnons in BiFeOs. 

PACS numbers: 75.40. Cx,75.85.+t,75.25.-j 
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INTRODUCTION 



The coupling between ferroelectric (FE) polarization and magnetic order in solid state 
materials, known as multiferroics, attracts significant interest because of its importance 
for understanding the properties of these materials, as well as of engineering application 
possibilities[l-4|. Among the known multiferroics, BiFeOs (BFO) is the only one exhibiting 
both magnetic and FE orders at room temperature. These orders are coupled, and control 
of both FE and magnetic properties by an electric field has been demonstrated in thin 
film and single crystal samples, as well as discussed theoretically. [7|, |8| Importantly, 
the magnitude of the spontaneous electric polarization (P ~ lO^/xC/cm^) in BFO is very 
large, facilitating control by the electric field, and opening possibilities for applications 
in multistate memory devices, spintronics, and configurable electronics. 9|] Therefore, BFO 
attracts considerable attention compared to other multiferroics. Until recently, BFO research 
was limited to thin films showing rather different properties from those of the bulk samples, 
and to polycrystals. Since suitable single crystals have become available recently, their 
properties were subject of many studies g, ^1^. Some of the basic properties of BFO, 
however, still remain to be characterized. 

The temperature-dependent evolution of the magnetic structure in BFO is an important 
subject with several currently unresolved issues. Below the FE transition at T^. ~1150 
K, BFO exhibits an R3c rhombohedral perovskite structure, which is described in this 
paper using a pseudo-cubic notation with a ~3.96 A, and a ~89.4°. Fe^"*" spins order at 
T/v ~640 K. It has long been known that the magnetic order is of the antiferromagnetic 
G type, with a long-range cycloidal modulation superimposed, [isl This modulation can 
propagate along three directions equivalent by symmetry, described by propagation vectors 
n = (5(1,-1,0), T2 = 5(1,0,-1), and rg = 5(0,-1,1), where 6 ~0.0045 r.l.u. For each of the 
three cycloids, the spins rotate in the plane defined by the (111) and r vectors, see Fig. 
l(a,b). A possibility for a so far undetected canting of the spins out of the cycloid plane has 
been discussed theoreti cally , j? | Below the room temperature, several anomalies have been 



observed in Raman 
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161 ] and dielectric susceptibility measurements. [17] Some of these 
anomalies, such as those at T=140 K and 200 K, are quite significant. It was proposed 
that these anomalies were associated with putative spin reorientation transitions, which 
are rather common in magnetic orthoferrites.|l| We note that spin reorientation can occur 
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through magnetic transitions within the existing magnetic domains, or through changes 
in the fractional populations of the equivalent magnetic domains Tj. NMR measurements 
indicate that the spin cycloid is anharmonic at T=5 K, while the anharmonicity is absent 



at T=300 K. 
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g|-|20| However, this effect is subtle, and therefore unlikely to produce the 
strong anomalies at T=140 K and 200 K. In any case, the NMR measurements of the 
magnetic order are indirect, and require confirmation. In fact, neutron diffraction studies 
of polycrystals did not find any measurable differences between the magnetic structures at 
T=4 K and at the room temperature. 21j In particular, no dependence of 6 on temperature 
was reported. So far, measurements of the magnetic order parameter and search for possible 
changes of the magnetic structure by direct techniques in the full temperature range, from 
helium temperatures up to T^, were absent, in part due to absence of suitable single crystals. 
Such measurements are clearly necessary to address the questions discussed above. They are 
also needed to understand the character of the magnetic transition at T/v and its possible 
interrelationship with the FE properties of BFO. The neutron diffraction and magnetization 
measurements described in the current study were motivated by these needs. 



EXPERIMENTAL DETAILS 

Single crystals of BiFeOs were grown using flux method as described in Refs. 
The crystal growth occurs below the FE transition temperature and tends to produce FE 
monodomain crystals. Neutron diffraction experiments were carried out on BT-9 triple-axis 
spectrometer at NIST Center for Neutron Research. A closed cycle cryostat and a furnace 
were used for temperature-dependent measurements. The neutron energy was 14.7 meV, 
and collimations of 40-40-S-40-open (order parameter measurements), or 10-10-S-lO-open 
(third harmonic studies) were used. To reduce the A/2 contamination of the neutron beam, 
pyrolytic graphite filters were used before and after the sample, as well as in-pile. A 3x2x0.5 
mm^ FE monodomain sample was oriented in the [h, k, {h + k)/2) scattering plane. In this 
paper, we use the pseudo-cubic notation for all the reciprocal lattice vectors, which are 
measured in reciprocal lattice units, r.l.u. Magnetization measurements were performed 
using a SQUID magnetometer under an applied magnetic field of 0.2 Tesla. In all figures, 
error bars are from counting statistics. 
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RESULTS AND DISCUSSION 



To characterize temperature-dependent behavior of the magnetic order parameter, pe- 
riodicity of the spin cycloid, equivalent magnetic domain populations, and possible anhar- 
monicity of the cycloid, single-crystal neutron diffraction experiments were carried out. The 
most intense magnetic reflection is observed at the Qq ± positions, where Qo=(0.5, 0.5, 
0.5), and i=l,2,3, giving rise to six diffraction peaks. The experimental scattering geometry 
is illustrated in Fig. 1 (c,d). The horizontal scattering plane contains perpendicular vectors 
Qo and Ti. Therefore, magnetic domain ti (the in- plane domain) produces peaks at Qo±S 
in the scan along the Ti centered at Qq. The magnetic domains T2 and T3 (the out-of- plane 
domains) produce peaks in the plane normal to the scattering plane. Nevertheless, they 
are also captured by this scan because of the broad instrumental resolution normal to the 
scattering plane. As illustrated in Fig. 1(d), the out-of-plane domains produce overlapping 
peaks at Qq ±S/2. Thus, a single scan along ri at the (0.5, 0.5, 0.5) position allows to char- 
acterize the in-plane domain, as well as to measure the combined signal of the out-of-plane 
domains. 

The scans were taken in the temperature range from T—5 K to 650 K. The scattering 
profiles are shown in Fig. 2 for several representative temperatures. For T > Tn ~640 
K, a weak temperature-independent peak is observed at Qq. It can originate from a A/2 
leakage in the neutron beam, or reflect a previously undetected small lattice distortion at 
this wave vector. In the magnetic phase, two asymmetric peaks centered roughly at Qq±S 
are observed. These peaks have clear shoulders in the direction of Qo, as best seen for the 
T—5 K data in Fig. 2. The overall peak shape is described perfectly by combined signals 
from the in- and out-of-plane magnetic domains at Qo^S, and Qo±6/2, respectively. This 
picture was confirmed by the room-temperature scans in the planes containing Qq and either 
T2 or Ta (not shown). Therefore, the scattering profiles were fitted using five Gaussian peaks: 
at Qo^S (the in-plane magnetic domain), Qq±S/2 (combined out-of-plane domains), and a 
weak temperature- independent background peak at Qq. The common Gaussian peak width 
was determined as a global parameter fitted to a common value over the entire temperature 
range. The intensities of the magnetic peaks, as well as the wave vector of the cycloid 5, 
were fitted at each temperature. This procedure resulted in high-quality fits, as shown in 
Fig. 2. 
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Figure 3(a) shows the temperature dependence of the integrated magnetic intensity, which 
is proportional to the square of the magnetic order parameter. It exhibits a canonical 
second-order-like behavior, and is well described by the power law in a rather 

large region below Tjy, see Fig. 3(a). The power law gives T/v=639.6(5) K, and /3=0.34(3). 
The critical exponent (3 is in agreement with three-dimensional magnetic transition. A 
similar /3 value was previously found in Mossbauer studies. 22|] The transition temperature 
Ttv is within the range of previously reported values for polycrystalline samples (590-650 
K),[l|, I23I and is consistent with theoretical predictions. 24] The ratio of the integrated 
intensity of the magnetic peaks at Qq±6 to that at Qq±6/2 gives the ratio of the in-plane 
and the out-of-plane domain populations. This ratio is shown in Fig. 3(b) as a function of 
temperature. Except for the narrow vicinity of the magnetic transition temperature, this 
ratio is essentially constant. As T/v is approached, the ratio decreases from about 1.5 to 
0.8, i.e. the relative population of the out-of-plane domains grows. This observation may 
reflect increased significance of thermal fluctuations for T ~ T/y. The magnitude of the 
cycloidal modulation wave vector 6 decreases with temperature, as shown in Fig. 3(c). The 
corresponding modulation period grows from 629(5) A at T=5 K to 780(30) A at T=615 
K. There is an apparent change of slope of this temperature dependence at T ~400 K. In 
previous works, [21]] the measurements were done in a smaller temperature range, and no 
temperature dependence of 6 was reported. 

Figure 3(d) shows the temperature dependence of the magnetization, with the magnetic 
field applied along different crystallographic directions. The general behavior below T/v 
is consistent with an antiferromagnet containing a small ferromagnetic impurity fraction. 
Similar behavior was reported elsewhere. [l|, |25| The anomaly at T/v itself is, however, rather 
unusual. The abrupt upward magnetization jump which accompanies appearance of the 
magnetic order is seemingly inconsistent both with the otherwise second-order-like behavior 
of the magnetic order parameter, and with the antiferromagnetic transition. The magnetiza- 
tion jump is rather small, corresponding to effective magnetic moment of ~ 5 x 10~^yUB/spin 
in the applied field of 0.2 Tesla. This suggests that the spin component involved in the 
abrupt transition is small, which is consistent with the order parameter behavior of Fig. 
3(a). Interestingly, the distribution of the magnetic domain populations also exhibits an 
abrupt change in the same temperature region. These data suggest that the magnetic struc- 
ture of BFO is more complex than currently believed, or that nontrivial effects related to the 
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magnetic domains occur near the magnetic transition temperature. The symmetry of the 
magnetic structure is a key factor defining possible magnetoelectric effects, such as existence 
of hnear magnetoelectric coupling. [26j] Further study of the nature of the transition at T/v 
is, therefore, of significant interest. 

Deviations of the magnetic structure from the circular cycloid have been detected pre- 
viously, but only at cryogenic temperatures. Specifically, line shapes observed in NMR 
spin-echo measurements at T=5 K are inconsistent with this structure. 18|-|20[| The NMR 
results were explained using an anharmonic cycloidal spin structure, which results from 
magnetic anisotropy in the plane of the spin rotation. Free energy density in BFO can be 
written [26| in the form / = fh + fexch + fan- Here is the relativistic inhomogeneous magne- 
toelectric interaction, fexch = ^Y^iSUY (i = x,y, z) is the inhomogeneous exchange energy 
with exchange stiffness A, fan = KuSin'^{9) is the anisotropy energy with the anisotropy 
constant K^-, 1 is the antiferromagnetic vector, and 9 is the angle between 1 and the (111) 
direction. The solution of this model is an anharmonic cycloid in which the spin component 
along the (111) direction is given by the elliptic Jacobi function sn{AK{m)x / \,m), where 
m is its parameter, K{m) is the complete elliptic integral of the first kind, A is the modu- 
lation period, and x designates the propagation direction. [26] Depending on the parameter 
m, the spin modulation changes from purely circular (m=0), to a square wave in which the 
spins bunch along the (111) direction (m=l). NMR line shape analysis based on this model 
indicated that the spin cycloid is strongly distorted at T=5 K (m=0.95), while it becomes 
essentially circular at room temperature (m=0.48). 

Magnetic neutron diffraction signal is produced by the spin components normal to the 
scattering vector. In our experiment, the scattering vector is essentially parallel to (111), 
and therefore only the spin components along the cycloid propagation vector r contribute 
to scattering. These components are given by (l-sn^(4i^(m)x/A, m))°-^=cn(4ii'(m)x/A, m). 
The diffraction intensity is proportional to the square of the Fourier transform of the mag- 
netization density. The Jacobi function has only odd harmonics in its Fourier series, and 
therefore scattering from the anharmonic cycloid described above should produce peaks at 
Qo ± nr {n is odd), whose intensity is proportional to the square of the appropriate Fourier 
coefficient. The ratio of these coefficients depends only on m. Thus, the ratio of the in- 
tegrated intensities of the first (^=1) to the third {n=3) harmonics of the magnetic peak, 
/1//3, provides direct measurement of m. Fig. 4(a) shows high-resolution scans capturing 
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the Qo — Ti and Qq — 3ri positions, taken at T=5 K and 300 K. A barely discernible third 
harmonic appears to be present at T=5 K, while none is seen at T=300 K. The data in 
Fig. 4 (a) were fitted using the model described above with the third harmonic peak added, 
and the line shapes changed to the resolution-corrected pseudo-Voigt function to correctly 
describe the wings of the peaks. For T=5 K, the fit gives /i//3=500±100 corresponding to 
m=0.50, and for T=300 K /1//3 >2000. The obtained m=0.50 produces only a very small 
deviation from the circular cycloid structure, as illustrated in Fig. 4(b), and in the inset in 
Fig. 4(a). In contrast, the value of m=0.95 derived from the NMR data at T=5 K would 
imply a clear bunching of spins along the (111) direction, see Fig. 4(c). We note, however, 
that m=0.95 would give /1//3 ~25, which is clearly inconsistent with our results. Neutron 
diffraction provides a more direct method for measuring the anharmonicity than NMR, and 
therefore we conclude that the low-temperature cycloid is more harmonic than previously 
reported. Both of these techniques indicate that the cycloid is essentially circular at room 
temperature. 

The anisotropy energy Ku is related 19[ to m by Ku = 16mAK'^('m)/X'^. The literature 
values of the exchange stiffness A are in the range of (2-4) x 10"^^ J/ cm. 18|, ll9|, l27H30| Using 
m=0.50, one obtains -ft'„=(1.5-3) xlO"^ J/cm^=(0.6-1.2) xlO"^ eV/spin. This is about four 
times smaller than that obtained using the NMR value of m, and two times smaller than 



the value found using ESR. 



3l| We note that the magnetic anisotropy is comparable to 



the inhomogeneous exchange energy E^xch ~ A6'^ since Ku/{A6'^) = AmK'^ (m) / tt'^ =0.7 for 
m=0.50. Both of these energies are smaller than kT, and therefore the anharmonic cycloid 
can be significantly affected by thermal fluctuations. This could explain, in particular, 
the apparent disappearance of the third harmonic peak at high temperatures, since the 
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fluctuation-induced peak broadening often increases with the harmonic order, 
any broadening of the already small peak would render it immeasurable. 

Anharmonicity of the spin cycloid is an important parameter defining dynamical magneto- 
electric response. It determines which spin wave branches couple to the electric polarization, 
giving rise to hybrid excitations (electromagnons) . It was predicted, for instance, that the 
cyclone magnon mode 02? and the extra-cyclone magnon mode ijj^ only have the hybrid 
character in the presence of the third harmonic in the cycloid. 3J] Higher-order harmonics 



make other cycloidal spin modes hybrid. Such spin modes have been observed in BFO. 
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Importantly, they have been found to exhibit giant energy shifts in an applied electric field, 
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indicating a significant magnetoelectric coupling. |14j| Our measurements of tlie cycloid an- 
harmonicity and of the relevant magnetic interaction energies should provide valuable input 
for microscopic description of the excitation modes and of the observed giant magnetoelectric 
effects in BFO. 

As discussed above, various anomalies have been observed below room temperature in 
Raman and dielectric property studies. Some of these anomalies were attributed to putative 
magnetic transitions. We do not observe any changes in the magnetic structure for T <300 
K, except for a slow and gradual change of the cycloid periodicity. Other recent neutron data 
are in agreement with these results. jsSi] Our data therefore show that the reported anomalies 
are not related to any magnetic transition, with only possible exception being a sudden 
change in the magnetic domain populations at these temperatures. While the temperature- 
dependent behavior of the domain populations could certainly be sample dependent, the 
observed absence of such a change in our samples makes the latter scenario appear unlikely. 
Other physical mechanisms should therefore be sought to explain the reported anomalies. 



CONCLUSIONS 



To summarize, we report neutron diffraction and magnetization studies of the magnetic 
properties of single- crystal BiFeOa in the temperature range from T=5 K to the temper- 
atures exceeding T/v. The magnetic order parameter exhibits a smooth second-order- like 
dependence typical of the three-dimensional magnetic transition. On the other hand, the 
magnetization exhibits a small jump at T/v, which is seemingly inconsistent with a regular 
antiferromagnetic transition. The reason for this behavior is yet to be determined. The 
period of the magnetic modulation smoothly increases with temperature from 630 A at 
T=5 K to ~780 A near T/v, exhibiting a change of slope at T ~400 K. The populations of 
the magnetic domains do not change with temperature, except in the close vicinity of T/v. 
There is a very small anharmonicity in the magnetic cycloid at T=5 K (the Jacobi function 
parameter m=0.50), which becomes undetectable at room temperature. It is significantly 



smaller than that reported in previous NMR studies. 



18H20| The magnetic anisotropy en- 



ergy derived from the observed anharmonicity is Ku ~10 eV/spin, it is comparable to the 



inhomogeneous exchange energy. The observec 
understanding the nature of hybrid excitations 



anharmonic parameters bear relevance for 



34l | as well as provide input for understand- 
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ing the giant effects [14!] of an electric field on the spin dynamics of BiFeOs. No magnetic 
anomalies have been found below the room temperature. We therefore conclude that the 
anomalies previously observed in Raman [l[ Q, 16] and dielectric measurements [ij] are not 



related to spin reorientation transitions. 
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FIG. 1: (a) Cycloidal magnetic structure of BiFeOs projected on the plane of the (111) and r 
vectors. About one tenth of the cycloid period along r is shown, (b) View of the same Fe chains, 
labeled 1 through 4, projected on the plane normal to (111), (c) Pseudo-cubic unit cell of BFO 
and the magnetic modulation wave vectors n originating at Qo={0.5, 0.5, 0.5). (d) Scattering 
geometry. The horizontal scattering plane contains ri and (111). Circles indicate the magnetic 
peak positions. The experimental resolution function is sketched with a dark vertical slab. 
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FIG. 2: Magnetic diffraction peaks in the vicinity of Qo={0.5, 0.5, 0.5). The scans are centered at 
Qo, and the scan direction is along (H,-H,0.5). Solid lines are results of fits discussed in the text. 
Dashed and dotted lines depict the calculated contributions due to the in-plane and out-of-plane 
magnetic domains, respectively. 
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FIG. 3: Temperature dependences of (a) Magnetic order parameter given by the integrated scatter- 
ing intensity near Qq, (b) Ratio of the population of the in-plane magnetic domain to the combined 
populations of the out-of-plane domains, (c) The magnitude 5 of the wave vector (5(1,-1,0) of the 
magnetic modulation, and (d) The magnetization in the field of H=0.2 Tesla for various directions 
of H (1 emu/Oe = 47rxl0^^ m^). The solid line in (a) is the result of the power-law fit discussed 
in the text. 
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FIG. 4: (a) Magnetic diffraction peak profiles for T=5 K and 300 K on a logarithmic scale. Solid 
lines depict results of fits discussed in the text. Dashed lines indicate positions of the first and third 
harmonics. The r=300 K scan is shifted down for clarity. The inset shows the spatial variation of 
the spin component Sr along the modulation propagation direction r for the circular cycloid with 
m=0 (top curve), and for the anharmonic cycloid at T=5 K with m-=0.50 (bottom curve), (b) Fe 
spins in the anharmonic cycloid with experimentally determined m=0.50 plotted from the same 
origin. For clarity, only one third of all spins are shown, (c) The same as (b), but for m=0.95. 
Our results are inconsistent with the latter value. 
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